Electron- and phonon transport in silicon nanowires: an atomistic approach to 

thermoelectric properties 
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We compute both electron- and phonon transmissions in thin disordered siUcon nanowires. Our 
atomistic approach is based on tight-binding and empirical potential descriptions of the electronic 
and phononic systems, respectively. Surface disorder is modeled by including surface silicon vacan- 
cies. It is shown that the average phonon- and electron transmissions through long SiNWs containing 
many vacancies can be accurately estimated from the scattering properties of the isolated vacancies 
using a recently proposed averaging method [Phys. Rev. Lett. 99, 076803 (2007)]. We apply this 
averaging method to surface disordered SiNWs in the diameter range 1 — 3 nm to compute the 
thermoelectric figure of merit, ZT. It is found that the phonon transmission is affected more by 
the vacancies than the electronic transmission leading to an increased thermoelectric performance 
of disordered wires, in qualitative agreement with recent experiments. The largest ZT > 3 is found 
in strongly disordered (111) oriented wires with a diameter of 2 nm. 

PACS numbers; 63.22.Gh, 66.70.-f, 73.63.-b, 72.10.Fk 



I. INTRODUCTION 

Both electronic and thermal transport in quasi one- 
dimensional systems such as carbon nanotubes and semi- 
conducting nanowires are of great importance for their 
apphcations. Often a large electronic conductance, Ge, 
and a large thermal conductance, Gph is desired in elec- 
tronic components. However, for thermoelectric applica- 
tions, the thermal conductance should be as small as pos- 
sible to increase the efficiency. The performance of a ther- 
moelectric device is characterized by the figure of merit, 
ZT = S'^ctT/k, where S is the Scebeck coefficient, a the 
electronic conductivity, T the temperature, and k is the 
total thermal conductivity. A good thermoelectric device 
should have ZT > 3 The ideal thermoelectric mate- 
rial is a so-called 'phonon-glass-electron-crystal' (PGEC) 
with a high thcrmopower, S*^ a and a low thermal con- 
ductivity, K. Different approaches have been followed 
to increase ZT including supcrlatticcs, alloys, embedded 
nano-particles in both 3D, 2D, and ID structures 0. It 
was shown theoretically in Ref. 3 that the thermoelec- 
tric figure of merit is larger in one-dimensional structures 
than in bulk. Indeed, recent experiments showed very 
large figure of merits in silicon nanowires (SiNWs) 
even though bulk silicon is a poor thermoelectric mate- 
rial. These experiments indicated that the phonon con- 
ductance was reduced more by surface disorder in SiNWs 
than the electronic conductance, leading to very high ZT. 

From a theoretical point of view it is not obvious that 
surface roughness will affect the phonons more than the 
electrons, and a model that treats both effects on the 
same level of approximation is thus desired. In order 
to theoretically interpret experiments or possibly to de- 
sign new thermoelectric materials or devices, it is thus 
important to treat both the electronic and the phononic 
transport on the same footing. 



Recently, Vo et al. [Gj] used ab initio calculations and 
the Boltzmann equation to calculate ZT, and found val- 
ues up to 8 in SiNW. The phonon heat conductance was, 
however, not calculated but used as a free parameter. 
However, we recently showed that the phonon heat con- 
ductance of pristine SiNWs is strongly anisotropic with 
(110) wires having up to two times larger conductance 
than (100) and (111) wires 0]. We show in this arti- 
cle that the anisotropics remain in disordered wires, thus 
clearly effecting the thermoelectric properties of SiNWs. 

Phonon transport was recently modeled in nanowires 
with surface disorder [1] and in SiNWs with an amor- 
phous coating 0. Also, several recent theoretical works 
have been concerned with scattering of electrons in 
SiNWs by surface rou g hne ss [lOl . [ll| and defects or 
dopant impurities [H, [H, O, Il5l|. We showed in Ref. [ll| 
that the average electronic conductance of an ensemble of 
long SiNWs containing many randomly placed dopants 
could be accurately estimated from the scattering prop- 
erties of the isolated dopants. This enables much faster 
calculations and the ability to study a larger diameter 
range. 

In this paper we demonstrate that the same single- 
defect averaging can be applied to both phonon- and elec- 
tron transport in the case of surface vacancy scattering. 
The average thermal conductance of long SiNWs con- 
taining many vacancies can thus be accurately estimated 
from the scattering properties of the isolated vacancies. 
We use this knowledge to compute both the electronic- 
and phononic conductances in ultra-thin SiNWs with sur- 
face disorder and to calculate the thermoelectric figure 
of merit in SiNWs with diameters ranging from 1-3 nm. 
We apply an empirical potential model to describe the 
phonons while the electronic system is modeled by a near- 
est neighbour tight-binding Hamiltonian. Both electron- 
and phonon conductances are calculated within the non- 
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equilibrium Green's function (NEGF) formalism [l6| . 

Our microscopic theory confirms the experimental 
trend [1]: the electrons are less affected by surface dis- 
order than the phonons, and the thermoelectric perfor- 
mance increases for increasing disorder (increasing num- 
ber of vacancies). We consider wires in the (100), (110), 
and (111) directions and find that (111) oriented wires 
have the highest thermoelectric figure of merit with an 
optimum diameter at 2.0 nm. To simplify our analysis, 
we neglect both phonon-phonon, electron-phonon and 
electron-electron scattering mechanisms. Due to the rel- 
ative simplicity of our model, our results do not allow 
for a quantitative comparison with experiments. How- 
ever, we believe that our results do shed light on both 
electron- and phonon surface scattering, and on the rel- 
ative magnitudes of the two. 

The paper is organized as follows: In Sect. |TT] we 
present the phonon empirical potential model and the 
electronic tight-binding model. Also, we explain how 
to calculate the electronic- and phononic transmissions 
through long wires either by a recursive method or, sim- 
ply from the transmission through single, isolated vacan- 
cies. In Sect, mil we present the results and end up by 
discussion and conclusions in Sect. IIVI 



unit cells, respectively. The calculations use periodic 
boundary conditions and the super cell method. Neigh- 
bouring, parallel wires are sufficiently separated and 
therefore do not interact. A pristine, defect free wire can 
be modeled using only a single unit cell in the super cell. 
We shall denote by koo and koi the dynamical matrices 
describing the pristine wire unit cell and the coupling be- 
tween neighbouring unit cells, respectively. Below, Fig. 
[3] illustrates the unit cells. 

When modeling vacancies, larger super cells arc needed 
in order to avoid interactions between the periodically re- 
peated vacancies. We find that 5 unit cells in the super 
cell H^l is sufficient for the scattering properties to be 
converged (see also Fig. [1]). We will denote the dynami- 
cal matrix around a given vacancy at position v by K''. 
When calculating the phonon properties we model pure 
silicon wires without any surface passivation. We have 
recently shown that including hydrogen on the surface 
only leads to insignificant changes Q. Also, we showed in 
Ref. that the thermal conductance of pristine SiNWs 
obtained with the TEP model agreed quantitatively with 
the more elaborate density functional theory (DFT) cal- 
culations. 



II. METHOD 

As a model of surface disorder we introduce surface sil- 
icon vacancies, i.e. removed surface silicon atoms. This is 
a very simple model and deviates from the standard de- 
scription of surface roughness, where the thickness of the 
wire diameter fluctuates on some characteristic length 
scale along the wire [ll|. However, the simplicity of the 
vacancy disorder allows us to study single vacancy scat- 
tering and apply the averaging methods of Ref. [13 and 
in this way study a larger diameter range than would be 
computationally manageable with our present methods 
and implementations. 

A. Phonon empirical potential model 

The phonon system is described by an atomistic model 
with the interatomic potentials parameterized by the Ter- 
soff empirical potential (TEP) model [l3| as implemented 
in the " General Lattice Utihty Program" (gulp) [H] . We 
use GULP to relax the atomic structure and to output the 
dynamical matrix, K, for the relaxed system. Since the 
Tersoff potential is limited to nearest neighbour interac- 
tions, K can be written in a block-tridiagonal form: 

/kn ki2 \ 

_ I k21 k22 k23 1 , . 

k32 k33 k34 ' 
V k43 k44/ 

where the sub-matrices k^^ and ki.i±i describe the force 
constants within a unit cell and between neighbouring 



B. Electronic tight-binding model 

The electronic system is described by a nearest neigh- 
bour sp^cf's* tight-binding (TB) basis. Contrary to the 
phonon system, it is necessary to include surface passiva- 
tion in the electronic description in order to saturate the 
Si dangling bonds, and to avoid localized states inside the 
band gap. Often hydrogen is used for passivation in DFT 
based calculations and we adopt this strategy by includ- 
ing a hydrogen-like atom with a single s-orbital. The 
purpose of the hydrogen is only to passivatc the dangling 
bonds, and not to describe the Si-H interactions in detail. 
The Si-Si TB parameters are from R ef. [2l]| while the Si- 
H interactions are the same as in Ref. [22l . The same Si-Si 
tight-binding parameters were recently applied to model 
surface roughness in SiNWs transistors [2J]. Due to the 
nearest neighbour model, the electronic Hamiltonian, H, 
has a block tridiagonal form as in Eq. ([T]) , in analogy with 
the dynamical matrix, K. When a Si surface vacancy is 
introduced, we passivate the resulting Si dangling bonds 
with extra hydrogen atoms. 



C. Electronic conductance 

We calculate the electronic current using the non- 
equilibrium Green's function (NEGF) formalism [23|. As 
usual in NEGF we divide our system into a left, central, 
and right region. The left and right contacts are modeled 
as two semi-infinite, defect free wires, and can be taken 
into account via the self-energies S£ j^(i?) and the 
retarded Green's function in the central region is calcu- 
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latcd as 

G^{E) = [EI-Uc- UliE) - T,},{E)]-\ (2) 

where He is the Hamiltonian in the central region con- 
taining defects, and I is the identity matrix. The trans- 
mission function through the central region is given by 

r{E) = Tr[TL{E) G-{E) Tn{E) G''(i?)], (3) 

wheieTL.R{E) ^(i;) - Sf, ^(i;)). The electronic 

current. / is given by the Landaucr formula, 

2p 

I=^J dETiE)[f{E,iiL)-f{E,^^R)], (4) 

where f{E, fiL.n) = 1/ (exp [{E - fiL,R)/kBT] + 1) is the 
Fcrmi-Dirac distribution function in the left and right 
leads with chemical potentials j-iL,R- The voltage dif- 
ference between left and right is V ~ {^l — HR)/e. The 
linear electronic conductance, Ge(/^) = dl / dV , is written 
as 

where we have introduced the function Lm {^J■) to be used 
later: 

D. Phononic thermal conductance 

The phonon transmission can be calculated in a mathe- 
matically similar way as the electronic transmission using 
the substitutions 

EI Lo^M (7) 
He ^ Kc (8) 

in Eq. Here M is a diagonal matrix with elements 
corresponding to the masses of the atoms and Kc is the 
dynamical matrix of the central region. The left and right 
contacts are again modeled as two semi-infinite wires 
with self-energies 'Sl^r{u!) [2^. Note that we use the 
same symbols for electron- and phonon Green's functions 
and sclf-cncrgics, but with different arguments: G{E) for 
electrons and G{u!) for phonons. We stress that using 
([5]) for is only valid when anharmonic (phonon-phonon) 
scattering can be neglected. The electronic analogue, 
Eq. ([3]), is likewise limited to mean- field theories [24| . 
In bulk Si, the room temperature anharmonic phonon- 
phonon relaxation length at the highest frequencies is 
^ai^max) ~ 20 nm and increases as Aa oc uj~'^ at lower 
frequencies Q. For relatively short disordered wires, the 
anharmonicity effect is thus of limited importance. 

The phonon thermal current at temperature T can be 
calculated from the transmission function as [2^, [2^, [27| 

Jph{T) = 2^ y AuwT{lo) [nB{TL) - nB[TR)], (9) 



where n_B(T) = (e'^'^/'^s^ _ is the Bose-Einstein dis- 
tribution and Tl,r = T ± AT/2. In the limit of smaU 
temperature difference AT, the phonon thermal conduc- 
tance Kph{T) = is 

(10) 

E. Modeling long wires by recursion 

Both the electron- and phonon transmissions through 
longer wires containing more vacancies at random posi- 
tions can be calculated using the recursive Green's func- 
tion method (RGF), as described in [3]. By repeatedly 
adding either pristine wire parts or parts from different 
vacancy calculations, a long wire is "grown". Ensemble 
averaged properties are found by repeating the calcula- 
tions for many wires with different configurations of the 
vacancy positions. Typically, several hundred calcula- 
tions are needed to converge the sample averaged prop- 
erties, which makes the long wire calculations very time 
consuming. 

A schematic illustration of a segment of a long wire 
with randomly placed vacancies is shown in Fig. [T] Each 
small box represents a wire unit cell and the three white 
circles represent three differently placed vacancies. As in- 
dicated by the shading, the two neighboring unit cells at 
each side of a vacancy arc slightly modified as compared 
to the pristine wire unit cell, represented by white boxes. 
The dynamical matrix (or Hamiltonian) inside each box 

||||||||||||N|||||||||||||| 

• • • c ■ • • • 

FIG. 1: Sketch of a wire segment containing three different 
vacancies, represented by white circles. Each box represents a 
wire unit cell, and the two neighboring unit cells at each side 
of a vacancy are slightly modified as compared to the pristine 
wire unit cell, represented by white boxes. 

is one of the diagonal sub-matrices in (U]), k^i and the 
coupling between two neighbouring boxes is described 
by the off-diagonal matrices, k^.i+i. We emphasize that 
when we model wires with more than one vacancy, we as- 
sume that the individual vacancies are separated in space 
and do not interact. This limits the minimum vacancy- 
vacancy distance to 2 nm in the axial direction. This 
method has recently been applied to calculate the elec- 
tronic transmissions through long SiNWs [3, [H, HI] . 

F. Modeling long wires from single defects 

Recently wc showed that the length dependent sample- 
averaged electronic transmission, T, of SiNWs with ran- 
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domly placed dopant atoms ('defects') could be accu- 
rately estimated from the transmissions through the sin- 
gle defects as [Uliil: 



A^dct + (1 - iVdcf) 



121 

To 



(11) 



where Ndcf = nL is the total number of defects in a 
wire of length L with defect density n. Tq is the ideal 
ballistic transmission of the pristine wire, and (T) = 
iJ2t=i is the average transmission of the M dif- 

ferent, isolated vacancies. Note that all the transmis- 
sions are energy (frequency) dependent. In the short wire 
limit, A^def — > the transmission equals that of a pristine 
wire T = %. For iVdef = 1, T = (T), while in the diffu- 
sive regime iVdef 3> 1, T cx The use of Eq. pl|) is 
limited to the quasi-ballistic and diffusive regimes, where 
the inverse transmission, 1/T increases linearly with wire 
length. 
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FIG. 2: (color online). Thermal conductance and electronic 
transmission for 1.2 nm diameter (110) wires. Top row: Ther- 
mal conductance vs. length (a), at T = 300 K and vs. temper- 
ature (b) at L = 75 nm. Dashed line: Long wire calculations 
and sample averaging. Solid line: results using Eq. Bot- 
tom row: Electronic transmission vs. energy for holes (c) and 
electrons (d). Solid line: pristine wire transmission; single 
vacancy estimates with A^vac = 5 (dotted red) and A'^vac = 20 
(dashed blue). The markers show sample average results for 
A*'vac ~ 5 (red circles) and A'vac = 20 (blue squares). 

Figure [5] shows the thermal conductance vs. length 
(a) and vs. temperature (b). The results are obtained 
for a 1.2 nm diameter (110) wire with a mean vacancy- 
vacancy separation of 3.8 nm. The dashed line is ob- 
tained by sample-averaging the transmission of 200 wires, 
each 150 nm long (400 unit cells) with different vacancy 
positions, while the solid lines are obtained from Eq. 
PT|) . It is evident that the sample-averaged conductance 
can be accurately estimated from the single- vacancy cal- 



culations. In this specific case, the number of matrix 
inversions required for the sample average calculations is 
> 10"^ times larger than the single vacancy-estimates. 

In Fig. [2] (c) and (d) we show the sample averaged 
electronic transmission for wires with iVvac = 5 (red cir- 
cles) and TVvac = 20 (blue squares) number of vacancies. 
The dotted red and dashed blue lines are obtained using 
Eq. pT|) with iVvac = 5 and A^vac = 20, respectively (sub- 
stitute iVdof with iVvac in Eq. pT|) ). Again wc see that 
the simple and fast single- vacancy method, Eq. (jll[) . pro- 
vides a reasonable estimate of the time consuming sample 
averaging results. The rest of the results in this paper 
are obtained from the electron and phonon transmissions 
through single, isolated vacancies using pT|) to calculate 
the length dependent transmission. 



G. Thermoelectric figure of merit 

The efHciency of a thermoelectric material can be char- 
acterized by the dimensionless figure of merit, ZT, given 

by m 



S^GeT 



(12) 



where S is the Seebeck coefficient, Ge the electronic con- 
ductance, T the temperature, Kph and the phonon- 
and electron contributions to the thermal conductance, 
respectively. Gg and Kp^ are given by Eq. (O and (|10p . 
while S and can be calculated when the electronic 
transmission T{E) is known [sol. [Sllls^: 



Ke(/i) 



1 ^i(m) 
eT Lo{n) 

1 
T 



L2{H) - 



(13) 
(14) 



where Lm{^i) is given by Eq. ([5]). The phonon contri- 
bution to S caused by phonon drag, is ignored in the 
present work, as our focus is on surface scattering, and 
inclusion of phonon drag is significantly more demanding 
involving clcctron-phonon coupling. Recent experiments 
indicated that a large phonon drag could be responsible 
for large ZT in SiNWs and our calculated S might 
be too low. 



III. RESULTS 

We consider SiNWs oriented along the (100), (110), 
and (111) directions. Cross sectional- and side views of 
the different wires are shown in Fig. [3] (top and middle 
panel). In the bottom panel we show the corresponding 
electronic band structure in the vicinity of the band gap. 
The corresponding phonon band structures can are re- 
ported in Ref. Q. In the following we show results for 
vacancy scattering in different wires. The transmissions 
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calculated from Eq. (|TT|) are averaged over all inequiva- 
lent surface vacancy positions. 
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FIG. 4: (color online). Average phonon transmission from 
Eq. (|lip through D = 2.0 nm wires containing 1 (dotted red) 
and 10 (dashed blue) vacancies. The solid black curve shows 
the pristine wire transmission. 



FIG. 3: (color online). Cross sectional (top panel) and side 
views (middle panel) of (100), (110), and (111) oriented wires 
(from left to right). The wires are oriented along the z- 
axis. Note that we have not plotted the passivating hydrogen 
atoms. Bottom panel: Electronic band structure around the 
band gap. 



spond to hole- and electron transmissions in the valence- 
and conduction bands, respectively. Generally, the holes 
are scattered more than the electrons, and we would thus 
expect n-type wires to have higher ZT values than p-type 
wires. 



A. Phonon- and electron vacancy scattering 

Figure 5] shows the average phonon transmission func- 
tion through wires containing surface vacancies. The 
results arc calculated from the single vacancies using 
Eq. ini) with A'^dcf = 1 (dotted red) and iVdcf = 10 
(dashed blue). The solid black curve is the pristine wire 
{Ndci = 0) transmission. Observe that the four acoustic 
modes at fiio < 4meV are unaffected by the vacancies 
and transmit almost perfectly. Generally, the scattering 
is relatively weak in the interval 30 < < 40 meV, while 
a significantly stronger scattering is observed between 10- 
30 meV. Notice also that the (110) wire has the largest 
transmission while the (111) wire has the smallest. This 
anisotropy was recently analyzed in Ref. 0] for pristine 
wires. Based on the anisotropic phonon transmissions in 
Fig. [3] we would expect (111) wires to have the highest 
ZT values since they have the lowest thermal conduc- 
tance. On the other hand, (110) wires are expected to 
have the lowest ZT values since they have the highest 
thermal conductance. 

Figure [5] shows the average electronic transmission 
from Eq. (jlip . again with TVdof = 1 (dotted red) and 
-^dof = 10 (dashed blue). The solid black curve shows the 
pristine wire transmission. Left and right panels corre- 
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FIG. 5: (color online). Average hole- and electron transmis- 
sion from Eq. (|lip through D = 2.0 nm wires containing 1 
(dotted red) and 10 (dashed blue) vacancies. The solid black 
curve shows the pristine wire transmission. 
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B. ZT calculations 

Figures |4] and [5] indicate that both electrons and holes 
generally are less affected by the vacancies than the 
phonons, implying that increasing the surface disorder 
(increasing -/Vdet either by increasing the vacancy density 
or the wire length) would increase the ZT values. This 



N 



0.6 
0.4 
0.2 



(C) 




(b) 



(d) 



<100> 



(11 0> 



(f) 



<111> 



-0.6 -0.4 -0.2 1.3 1.6 1.9 2.2 

^ (eV) 



FIG. 6: (color online). Calculated ZT values vs. chemical 
potential at T = 300 K. Solid red: iVvac = 10; dashed blue: 
iVvac = 100; dotted black: iVvac = 1000. 

is indeed what we observe in Fig. [S] showing calculated 
ZT values as a function of chemical potential at different 
degrees of disorder. The temperature is T = 300 K. Both 
for electrons and holes in all wire directions, we observe 
an increase in ZT as the number of vacancies is increased 
from 10 (sohd red) to 100 (dashed blue). At TVvac = 1000 
(dotted black) the ZT is increased in all n-type wires 
(electrons) while only the (111) wire shows increased hole 
ZT . As expected, the (111) wires have the largest ZT 
values, while the (110) wires have the smallest. This is 
partly due to the anisotropic phonon heat conductance, 
and partly due to difference in the electronic band struc- 
ture. While both the (100) and (110) wires have the con- 
duction band minimum (CBM) at fc = 0, the (111) wire 
has an indirect band gap, with the CBM at fc ^ 0. This 
implies that close to the conduction band edge, there is 
a peak with six conducting channels in the (111) wire. 
The (100) wire quickly reaches a transmission plateau 
of four, whereas the (110) wire has a relatively broad 
plateau with only two channels. 

The maximum ZT occurs at a chemical potential close 
to the band edge, which requires a high doping concen- 
tration. The dopant atoms would certainly affect the 
electronic conductance [l^ [H, HI] and, depending on 
the dopant type, also the phonon transmissions. Us- 
ing phosphorous or aluminum as n- and p-type dopants 



would probably affect the phonons relatively little due 
to the small atomic mass difference between P/Al and 
Si, whereas nitrogen or boron are expected to scatter 
the phonons relatively more. As our main focus in on 
the surface disorder scattering the present calculations 
do not take the dopant scattering into account, and the 
computed values for Ge give an upper limit. 

Figure [7] shows the maximum ZT vs. wire diameter, 
D, for (111) wires with different degrees of disorder (num- 
ber of vacancies) . For both n-type (electrons) and p-type 
(holes) pristine wires with A^vac — 0, the maximum ZT 
increases with decreasing diameter. The reason is that 
at T = 300 K effectively all phonon modes contribute to 
the phonon heat conductance, which therefore scales as 
in the continuum limit, n^h oc A, where A is the cross 
sectional area. The electronic conductance, Ge, is more 
or less constant as long as the distance between the elec- 
tronic bands is larger than fc^T, which is the case for the 
diameter range considered. For larger wires, Ge will also 
increase. 

When vacancies are introduced, wc observe an opti- 
mal diameter at 2 nm for n-type wires (electrons) in- 
dependently of disorder strength For diameters > 
20 A, we also observe that increasing the disorder in- 
creases the maximum ZT, while for the smallest wires 
the largest ZT is found for A^vac = 100. At D = 20 A, 
and A'vac = 1000, the maximum ZT is six times larger 
than in the pristine wire with the same diameter. In 
the p-type wires (holes) with Avac = 10 and A^vac = 100 
the maximum ZT increases for decreasing wire diameter. 
For the strongest disorder (A^vac = 1000) we again see a 
maximum at Z? = 20 A, where ZT is 18 times larger than 
in the pristine wire. 
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FIG. 7: (color online). Diameter dependence of the maximum 
ZT in (111) wires with different numbers of vacancies 
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IV. DISCUSSION AND CONCLUSION 

It is interesting that an optimal diameter exists where 
ZT is maximized. To understand this we note that the 
weight of electron Bloch states, j^/'/d^, and phonon eigen- 
modes, |tip, on the surface atoms have a different diam- 
eter dependence as shown in Fig. [5] The weight of the 
lowest conduction band electronic Bloch states (in the 
F-point) on the surface atoms decays as ^ D^^ (circles). 
We note that this is the same scaling as one obtains from 
the effective mass Schrodinger equation. The average 
amplitude of all the phonon eigenmodes (also in the F- 
point) on the surface atoms scales, on the other hand, as 
~ (squares). 




16 25 
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FIG. 8: (color online). Diameter dependence of the average 
weight of the lowest conduction band electron Bloch state, 



\'^k{raurf)\ (circles), and average of all phonon eigenmodes, 
\u{rsurf)\^, (squares) 
tice the log-log scale. 



|w(?"su7-/)|^, (squares) at the surface atoms in (111) wires. No- 



In a Fermi's golden picture, the electron scattering rate 
7e oc \'4'k{i'vac)\'^ of & vacancy at position r^ac, and like- 
wise for the phonon scattering oc |u(ri;ac)P 01- The 
different diameter scalings of |V'fc(''sMr/)P and \u{rsurf)\^ 
imply that surface disorder affects the electrons relatively 
more for the thinner wires. This effect tends to de- 
crease ZT for decreasing diameter. On the other hand, 
as discussed above, the phonon conductance scales as 
oc thus increasing ZT for decreasing diameter. 
Given the two opposing diameter dependencies, it is rea- 
sonable that an optimal diameter exists where ZT is 
maximized. 



Summariziiig, we have used the averaging method pro- 
posed in Ref. [1^ and have shown that both the electron- 
and phonon transmission through long wires contain- 
ing many randomly placed vacancies can be reproduced 
from the transmissions through isolated vacancies. This 
enables one to perform relatively fast calculations on 
wires with experimentally relevant lengths. The aver- 
aging method was in Ref. (l5| applied to scattering at 
dopant impurities, and with the present study of vacan- 
cies, we believe that its validity is general for any random 
arrangement of localized structural disorder (defects, im- 
purities, vacancies, adatoms, physisorbed molecules, etc.) 
and for both electron- and phonon transmissions. It re- 
mains to be investigated whether a similar approach can 
be extented to systems where the wire diameter fluctu- 
ates along the wire. 

SiNWs oriented along the (111) direction have the 
largest ZT values while (110) wires have the smallest. 
This is primarily due to the anisotropic heat conductance 
m SiNWs 0. (110) wires have larger Uph and thus a 
smaller ZT . Differences among the wire orientations in 
the electronic transmission and Seebeck coefficients do 
also play a role. 

For the (111) wires, we observe an increasing maximum 
ZT for decreasing the diameter down to D = 2.0 nm 
which indicates that decreasing the wire diameter from 
typical dimensions of 10-50 nm to the sub-10 nm range 
might increase the thermoelectric performance. More 
thorough analysis including more realistic surface rough- 
ness is however required before firm conclusions can be 
drawn. We will address this in a future communication. 
Also, a quantitative model should also include electron- 
phonon- as well as phonon- phonon scattering [ssj . How- 
ever, we believe that the present results, where electron- 
and phonon-vacancy scattering is based on the same 
atomistic structure, is an important first step on the way. 
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